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Cosmological parameter estimation exercises usually make the approximation that the three stan- 
dard neutrinos have degenerate mass, which is at odds with recent terrestrial measurements of the 
difference in the square of neutrino masses. In this paper we examine whether the use of this approx- 
imation is justified for the cosmic microwave background (CMB) spectrum, matter power spectrum 
and the CMB lensing potential power spectrum. We find that, assuming Smls ~ 2.5 X 10"^eV^ m 
agreement with recent Earth based measurements of atmospheric neutrino oscillations, the correc- 
tion due to non-degeneracy is of the order of precision of present numerical codes and undetectable 
for the foreseeable future for the CMB and matter power spectra. An ambitious experiment that 
could reconstruct the lensing potential power spectrum to the cosmic variance limit up to ^ ~ 1000 
will have to take the effect into account in order to avoid biases. The degeneracies with other pa- 
rameters, however, will make the detection of the neutrino mass difference impossible. We also show 
that relaxing the bound on the neutrino mass difference will also increase the error-bar on the sum 
of neutrino masses by a factor of up to a few. For exotic models with significantly non-degenerate 
neutrinos the corrections due to non-degeneracy could become important for all the cosmological 
probes discussed here. 

PACS numbers: 



INTRODUCTION 

Standard cosmological measurements offer an excellent 
probe of neutrino physics il| . In contrast to Earth based 
measurements of neutrino oscillations that measure the 
difference in the square of neutrino mass eigen states, the 
cosmology is sensitive to the absolute mass of neutrinos. 
At the moment, cosmology is the only viable alternative 
to the beta-decay experiments in this field and exceeds 
it in accuracy. 

The most natural and accurate way to measure neu- 
trino masses with cosmology is via measurements of the 
power spectrum of fluctuations in the cosmic microwave 
background (CMB) and the power spectrum of mat- 
ter fluctuations, which is one of the basic products of 
the galaxy redshift surveys. In the future, the non- 
Gaussianity introduced by the lensing of the CMB fluc- 
tuations by the intervening structures between the last 
scattering surface and us will prove to be an important 
tool to constrain the power spectrum of matter fluctua- 
tions and consequently neutrino masses Q |^ . 

Many of the recent parameter estimation papers have 
pu t bounds on the sum of the neutrino masses 0, B IS 
ESmiillllllllllllllll, with the upper bound on 
the sum of neutrino masses of the order 0.5eV-leV. The 
three neutrino species were assumed to be degenerate in 
mass. 

However, Earth based measurements of neutrino oscil- 
lations suggest not only that neutrinos have mass, but 
also that these masses are not equal |3i EH- Neu- 
trino that is produced in a flavour eigenstate, which is 
some linear superposition of mass eigenstates, can lat- 
ter be measured as being of a different flavour, due 
to different mass eigenstates acquiring a different phase 



during propagation. In particular, oscillations of atmo- 
spheric neutrinos suggest a squared mass difference of 
lAm^gl = |m§ - mil =~ 2.5 x IQ-^eV^ HSHJ, while 
solar neutrino observations, together with results from 
the KamLAND reactor neutrino experim ent, point to- 
wards Ato?2 = mi - =~ 5 X IQ-'^eV^ [lllii. Note 
that in the former case, the sign of the difference is not 
know, while in the latter it is. Given these constraints, it 
is possible construct two hierarchies of masses. Assum- 
ing the lightest neutrino to be of a negligible mass, the 
neutrino masses can be either in a so-called normal hier- 
archy with masses around ^ 0,~ -^/Attz^, ~ •^/Amjj 
or in an inverted hierarchy, in which case masses are 
^ 0,^ -y/lAmlgl , ~ ■y/jAmfgl. This means that it is 
possible to rule out the inverted hierarchy simply by 
measuring the sum of neutrino masses and excluding the 
^ m > 2y^ Ato23 region. 

Another way of distinguishing between the two hierar- 
chies is to try to measure the neutrino mass difference di- 
rectly with cosmology. The above results imply that the 
three neutrino families are non-degenerate at the level of 
~ 5% if the sum of neutrino masses is 0.5eV. Decreasing 
the sum to 0.2eV, the two neutrinos are non-degenerate 
at the level of 25%. Hence, it is timely to investigate 
what are the biases introduced by assuming degenerate 
neutrino masses in cosmological probes of neutrinos and 
this is exactly what this paper is set to do. 

Significant amounts of work in this direction was done 
before in . This paper focuses on the ability of future 
Large Scale Structure (LSS) and CMB experiments to 
constrain neutrino masses by performing the Fisher ma- 
trix analysis. They were the first to perform the numeri- 
cal integration of the CMB and LSS power spectra by in- 
dependently integrating more than one neutrino species. 
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We extend their work in several aspects. Firstly, we also 
include the lensing potential reconstruction in the set of 
datasets used to constrain neutrino masses. This has 
been studied in antoher recent paper j^, which shows 
that it would be impossible to distinguish between nor- 
mal and inverted hierarchies using even a very optimistic 
future experiment. Nevertheless, the authors of that pa- 
per do not use the multi-neutrino code and do not discuss 
degeneracy assumption further. Secondly, rather than 
expanding around a fiducial model, we study the general 
parameter space (although we also provide Fisher ma- 
trix analysis in order to check the effect of degeneracies) . 
Finally, we provide a confirmation of their results by an 
independent implementation of a linear code with more 
than one neutrino species. An analysis of future sensi- 
tivities of high redshift galaxy surveys and CMB data to 
measure neutrino masses and number of neutrino species 
can be found in ji^ . 

Assuming that the standard physics of the early uni- 
verse applies, the energy density and mass of a given 
neutrino species are related through 



94.2eV 



(1) 



Interactions in the early universe before neutrino decou- 
pling ensure decoherence and so flavour physics does not 
enter cosmology. On the other hand, changing energy 
density (or mass) of a given neutrino species, affects 
the CMB and matter power spectra in two ways. First, 
it changes the redshift of the matter-radiation equality, 
thus affecting the position and height of the peaks in the 
CMB power spectrum and the maximum in the matter 
power spectrum. Second, it damps the power spectrum 
on small scales, because relativistic neutrinos in the early 
universe behave as radiation and effectively stream away 
from over-dense regions. The free streaming wave- vector 
is given by 



0.01— ^Mpc 
lev 



(2) 



and the matter power spectrum is damped in scales k > 
kfs- For neutrinos of a small mass the second effect is 
considerably more important. 

In general, the neutrinos become non-relativistic at a 
redshift of 



2xl03f^) 



(3) 



This means that neutrinos lighter that about 0.5eV 
will become non-relativistic after the recombination and 
thus have a very small effect on the CMB fluctuations 
power spectrum. The neutrino mass can, however, be 
inferred from the gravitational lensing of the CMB fluc- 
tuations. The quantity of interest here is the projected 
gravitational potential (see e.g. 0|) 



(n) = -2 / *(nr) 



rr-i. 



-dr. 



(4) 



where r are conformal distances, is the conformal dis- 
tance to the surface of the last scattering and the integra- 
tion is performed along our past light cone. The spher- 
ical power spectrum of Cf*^ can be recovered using 
a variety of methods 28, 29l l30l Isil Is^ and essentially 
contains information similar to that of the matter power 
spectrum. 



NON-DEGENERATE NEUTRINOS 

If one relaxes the degeneracy assumption, the neutri- 
nos of different masses become non-relativistic at differ- 
ent times and have different free-streaming lengths, re- 
sulting in small corrections to the various power spectra 
discussed above. Since the difference in the squares neu- 
trino masses Am^^ is over an order of magnitude larger 
than Am^2, we will, for the time being, assume the lat- 
ter is zero. We thus have two neutrinos of the same mass 
TOi = m2 and the third one of a different mass 7713. We 
parametrise the masses in terms of the sum of neutrino 
masses ^ TOj and the fraction a of the total mass in the 
third neutrino mass eigenstate, so that 



TO3 



(5) 



This particular parametrisation has been chosen, be- 
cause it allows for both extreme possibilities, namely that 
all of the mass is in the third neutrino (a = 1) or that 
third neutrino is massless (a = 0) . The value of a = 1/3 
corresponds to the degenerate case. 

We use a modified version of the CAMB linear solver 
[3^ that can evolve two families of neutrinos of different 
masses separately |34j . The accuracy_boost parameter 
was set to 2, which should result in an accuracy around 
0.1%. In addition, we used transf ersJiigh_precision 
option. We have checked that the results are the same 
regardless of whether a full hierarchy integration is per- 
formed or a switch to series in velocity weight once neu- 
trinos become non-relativistic is used. Other parameters 
were set to their nominal values for a ACDM universe. 
The CMB power spectra discussed here were lensed using 
the algorithm discussed in j^. 

In Figure n we plot the lensed CMB power spectrum, 
the matter power spectrum and the projected gravita- 
tional potential power spectrum for three models con- 
taining either three massless neutrinos or three massive 
neutrinos with ^ — 2eV and a = 1/3 and a = 1. 
The figures correspond to the standard flat ACDM cos- 
mology. The energy densities of baryonic and cold dark 
matter as well as curvature were kept fixed so that hot 
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FIG. 1: The lensed CMB power spectrum (top), the mat- 
ter power spectrum (middle) and the lousing potential power 
spectrum (bottom). Solid line correspond to the standard 
ACDM model. Other two models plotted have ^ rrii — 2eV 
and a = 1/3 (dotted) or a = 1 (dashed). See text for discus- 
sion. 

dark matter component grows at the expense of cosmo- 

logical constant. The sum of neutrino masses that large 
is not compatible with current observations and was cho- 
sen for illustrative purposes only. The effect on CMB is 
largely determined by the change in the scale factor of the 
matter-redshift equality (aeq), the change in which affects 
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FIG. 2: This figure shows the change in the for a cosmic 
variance limited experiment (to £ — 2000) if one wrongly as- 
sumes degenerate neutrinos. The contours are at Ax^ of 1, 
5, 25 and 125 (from a = 1/3 line outwards). See text for 
discussion of other features on the plot. 



the CMB peak positions (due to change in the sound hori- 
zon) and heights (via the early integrated Sachs- Wolfe 
effect). However, the correction due to non-degeneracy 
assumption is very small. The matter power spectrum 
and the lensing potential power spectrum show very sim- 
ilar trends. If one neutrino contains all the mass, it has a 
smaller free streaming Icingth and consequently damping 
does not extend to scales as large as in the case where 
neutrinos have degenerate mass. However, because the 
total neutrino energy density is the same, the overall ef- 
fect is the same for kfs- The latter also implies that, 
contrary to what might be naively expected, the precision 
measurements at very small scales (such as those probed 
by Lyman alpha forest) will not be sensitive probes of 
neutrino mass differences. 



RESULTS 

How big are the discussed effects for realistic neutrino 
masses and experiments? To answer this question we 
calculated power spectra for a grid of models with ^ rrii 
between and leV and a between and 1. In each 
case the power spectra were calculated using both the 
approximation that we have 3 degenerate neutrinos with 
a given "^i ^nd the correct distribution of neutrino 
masses, evolving the two neutrino species separately. 

For the CMB power spectrum, there exist a natural 
limit for the accuracy with which the power spectrum 
can ever be measured. The so called cosmic-variance is 
a result of a finite number of spherical harmonic modes 
on the sky and is (for £ > 50) excellently approximated 
by a Gaussian distribution with an error given by 
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FIG. 3: Same as Figure|21but for the lensing potential recon- 
struction that is cosmic variance up to ^ = 1000. 



o-cv(Q) = '^e\J 1 

For each error in our parameter space we calculated the 
change in induced by the correction stemming from 
non-degeneracy, i.e. 

Ax' = (^^^-^-^- - f)^ (7) 

where Cf^dcgcn is the theoretical prediction if degeneracy 
is assumed, for an experiment that is cosmic-variance 
limited up to ^ = 2000. Such measurements could be 
performed with a future experiment with a beam size of 
4', temperature sensitivity of At ~ luK and polarisation 
sensitivity of Ap = -\/2At ~ 1.4/xK [Sg- The results are 
plotted in the Figure [2| This figure is worthy some dis- 
cussion. The horizontal axis is the sum of the neutrino 
masses, while the vertical axis correspond to the a pa- 
rameter defined above. The dashed line at a = 1/3 shows 
the degeneracy condition (where Ax^ is equal to zero by 
construction) and the vertical dashed line corresponds to 
the tightest limit on the sum of neutrino mass found in 
literature j9|, that is 0.4eV. A few recent preprints find 
tighter limits using the latest cosmological data: 0.3eV 
in and O.lTeV in [s^- Solid thin hues are contours 
of constant Ax^ as stated in the caption. The two solid 
thick lines correspond to the value of a required to sat- 
isfy the ATO23 = 2.5 X 10~'^eV^ condition. The upper 
branch corresponds to the normal hierarchy, while the 
lower branch corresponds to the inverted hierarchy. 

This results shows the expected result that the mag- 
nitude of correction due to the non-degeneracy is a func- 
tion of both relative mass difference and the total sum 
of the masses. If one takes earth-based measurements of 
the mass square difference seriously, then it seems that 
the standard degenerate neutrinos assumption is a good 



one: the approximation is either saved by being a too 
small relative effect at larger sum of masses or being a 
too small absolute effect anyway, when the sum of masses 
is small. However, if one is not bound by the small mass 
square difference (when considering exotic neutrino mod- 
els, for example), then the effect can be quite large and 
detectable with high precision in the future CMB exper- 
iments. 

A much more interesting picture emerges if one looks 
at the reconstructed lensing potential. The reference ex- 
periment discussed above could, ideally, reconstruct the 
lensing potential to a cosmic variance limit up to ^ = 1000 
[sgI . Figure 01 shows contours analogous to that of the 
Figure 121 but for the lensing potential instead. One can 
see that for ^ O.leV , the non-degeneracy correc- 

tions can become important at high confidence. Our esti- 
mate is conservative since it assumes that information on 
lensing potential is not available beyond £ — 1000. A real 
experiment will still have sensitivity to recover the lensing 
potential at i > 1000, albeit with a sub cosmic variance 
precision. We also note, that linear transfer functions 
were used, but it is unlikely that non-linear corrections 
would significantly destroy the sensitivity. 

Next we turn to the matter power spectrum. The 
power spectrum is more difficult to consider as it is not 
clear what an idealised experiment can do in the pres- 
ence of non-linear biasing and complications arising from 
astrophysical considerations. Therefore we calculate the 
relative change of the slope of the linear power spectrum 
at three nominal values of A: = 0.005/i/Mpc, 0.01/i/Mpc 
and O.l/i/Mpc (where h is the reduced Hubble's con- 
stant). These are plotted in the Figure01 We see that the 
changes in the slope of the linear power spectrum are of 
the order of a ~ 0.1% and very likely undetectable in the 
foreseeable future. Again, we note, however, that if the 
bound on the mass square difference is abandoned, the 
non-degeneracy correction can be quite large and exceed 
the 10% mark in some cases. 



FISHER MATRIX ANALYSIS 

In the above sections we have shown that in some parts 
of the parameters space the effect of the difference in neu- 
trino masses can produce sizeable differences. This 
implies that neutrino mass difference must be taken into 
account in order to avoid biases in the data. However, 
it does not necessarily mean that neutrino mass differ- 
ence will be detectable due to possible degeneracies with 
other parameters. In order to check for that effect we 
perform a Fisher Matrix analysis. We use the following 
parametrisation of the model 

0i = (^LUb, Wcdm, h, T, Hs, A, ^ m, , (8) 
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parameter Normal Hierarchy Inverted Hierarchy Degenerate 1 Degenerate 2 



Y,mi/eV 0.055 0.105 0.055 0.105 

a 0.95 0.043 1/3 1/3 



(Fa 



-1/2 



-1/2 



N^ = {a~\)/{F~^f/ 



0.1 
4.9 

2.2 

2 0.26 



0.1 
2.9 

1.2 

0.2 



(FEmE,„)-'''VeV 0.003 
^EWl^'smSm)-'/' 18.8 

{F^L^mf'^h'^ 0.08 
iV. = EW(i^E^E™)'/' 0.6 

TABLE I: This table shows the results of 



0.003 0.003 0.003 

36.6 18.8 36.6 

0.06 0.013 0.04 

1.7 4.2 2.5 

matrix analysis. See text for discussion. 



where parameters have their usual meaning in the cos- 
mological context. The Fisher matrix is given by 

where XY is either TT, EE, TE (temperature, E polari- 
sation and cross power spectra) and LL (lensing potential 
cross-spectrum). Since the exact experimental parame- 
ters are not the focus of this work, we simply assumed 
that the CMB TT power spectrum is known to cosmic 
variance for € < 2000, while TE, EE and lensing spectra 
are known to cosmic variance for ^ < 1000. 

The interpretation of the Fisher matrix is straightfor- 
ward: (i^ri)~^/^ is the expected 1 — a error on the mea- 
surement of the i-th parameters, assuming all other pa- 
rameters to he fixed. The value of {F^^Y/'^ gives the ex- 
pected 1 — (7 error on the measurement of the i-th param- 
eter taking into account possible degeneracies with other 
parameters while the direction of these degeneracies are 
given by the eigenvectors of F. A Gaussian nature of the 
posterior is assumed throughout. 

We have performed the Fisher ma- 
trix analysis, using the nominal values for 
most parameters, {uj\y,uJcd-m,h,T,ns, A) = 
(0.02,0.12,0.7,0.11,1.0,2.3x10-9). For the values 
of the remaining two parameters we took two represen- 
tative points, one for the normal and one for the inverted 
hierarchy, which also satisfy the atmospheric neutrinos 
constraint. We have also performed the analysis by 
fixing a = 1/3 and excluding that parameter from the 
analysis. The results are shown in the Table |l| 

There are a few interesting conclusions to be drawn. 
Firstly, we find a reasonable agreement (given the crude 
approximation of experimental performance) with for 
the degenerate case, where we find that the marginalised 
error on the sum of neutrino masses will be of the order 



0.013 — 0.04 eV. Secondly, we note that the degeneracies 
with other parameters completely destroy the detection 
of the neutrino mass difference. From a modest few sigma 
detection, the error increases several-fold. Thirdly, as- 
suming degeneracy severely decreases the accuracy with 
which the sum of neutrino masses can be measured. In- 
deed, the analysis of the eigenvectors of the Fisher matrix 
show that the two strongest degeneracies involve ^ m, a 
and h in one eigenvector and Wdm, ot and h in the other. 
The latter is yet another face of the degeneracy discussed 
in m. 

CONCLUSIONS 

In this paper wc have examined the importance of the 
degeneracy assumption, which is often used in literature 
to constrain the sum of masses of neutrinos from cosmol- 
ogy. Since the combination of the measurements of neu- 
trino oscillations from the earth-based experiments with 
the upper limit on neutrino masses from the cosmological 
experiments imply that the neutrinos are non-degenerate 
at the level of at least several percent, it is not obvious 
that this assumption is a justified one. 

By comparing model predictions for a model with 
two different neutrino masses with that of an equiva- 
lent model with three degenerate neutrinos we were able 
to show that the degeneracy assumption is indeed valid 
for the constraints based on the CMB and galaxy power 
spectra. In the case of CMB, this is true even for the cos- 
mic variance limited experiment, while the matter power 
spectra are unlikely to reach the accuracy required to de- 
tect neutrino mass difference due to various astrophysical 
constraints such as scale-dependent biasing. Moreover, 
the size of the effect is at the level of numerical precision 
of the present generation linear codes. 
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FIG. 4: This figure shows the relative change in the matter 
power spectrum slope at fc = 0.005/i/Mpc (top), 0.01/i/Mpc 
(middle) and 0.1/i/Mpc (bottom). The thin solid and dot- 
ted lines correspond to 0.1% (dotted), 1% and 10% difference 
and increase in value from the thin horizontal dashed line 
outwards. 

The best hope for detecting neutrino mass difference 
with cosmology lies in the CMB lensing potential recon- 
struction. In an idealised future experiment such as that 
of |36l |. the neutrino mass difference could bias the re- 
sults of parameter estimation for ^rrii ^ O.leV, since it 
significantly affects the x^- However, by invoking Fisher 
matrix analysis we have shown that it would be impossi- 



ble to measure it due to degeneracy with other parame- 
ters. These degeneracies could be broken by constraining 
ujdui and hin a independent manner. We have also shown 
that once the degeneracy assumption is lifted, the error 
on the measurement of the sum of neutrino masses is sig- 
nificantly worsened; even for only two neutrino masses 
discussed here. 

If there is no information on the neutrino mass differ- 
ence, as it is usually the case in exotic scenarios, such as 
those involving extra sterile neutrinos, then the correc- 
tions are potentially significant for all cases. The slope of 
the matter power spectrum can change as much as 10% 
and the change could be detected with high significance 
even for the CMB power spectrum if ^ > 0.4eV. 
However, recent analysis have shown that the number 
and masses of neutrino species are significantly con- 
strained even with the present generation cosmological 
data and thus we will not explore this further. 
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